Fast quantum limited read-out of a superconducting qubit using a slow oscillator 



Goran Johansson, Lars Tornberg, CM. Wilson 
Microtechnology and Nanoscience, MC2, Chalmers, 96 Goteborg, Sweden 

We describe how to perform fast quantum limited read-out of a solid state qubit biased at its 
degeneracy point. The method is based on homodyne detection of the phase of a microwave signal 
reflected by a slow oscillator coupled to the qubit. Analyzing the whole quantum read-out process, 
we find that the detection is indeed quantum limited and that this limit may be reached even using 
a resonance circuit with a low quality factor, thus enabling the use of short measurement pulses. As 
an example, we discuss in detail the read-out of a Cooper-pair box capacitively coupled to a lumped 
element LC-oscillator. Furthermore, we give formulas for the backaction while not measuring, and 
discuss optimal parameters for a realistic design capable of fast (~ 50 ns) single-shot read-out. 



Superconducting qubits are strong contenders in the 
race to build a quantum computer [l|. Accordingly, a 
great deal of effort has gone into developing fast and ac- 
curate single-shot read-outs for these devices. 

Speed is one important characteristic for a read-out. 
For accurate qubit detection, the read-out has to be faster 
than the qubit relaxation time. Moreover, to implement 
quantum error correction, the read-out must be quicker 
than the decoherence time. 

Another figure of merit, which characterizes the back- 
action, is the quantum efficiency, r\. According to fun- 
damental principles of quantum measurement, the mea- 
surement time, t ms , needed to distinguish the two qubit 
states is always longer or equal to half the dephasing 
time, t v , after which the qubit has lost its quantum 
coherence^. Thus, the quantum efficiency has an up- 
per bound of unity, rj = t v /2t ms < 1. 

Recently, there has been increasing interest in disper- 
sive read-outs for superconducting qubits 0, 0, 0- 
which measure the reactive response of an oscillator cou- 
pled to the qubit. The energy needed for detection 
is then dissipated far away from the qubit, giving the 
schemes low backaction. Also, many dispersive read-outs 
work with the qubit biased at its degeneracy point where 
the decoherence induced by low frequency fluctuations is 
minimized 

In this letter, we consider how to optimize a dispersive 
qubit read-out for speed, while maintaining high quan- 
tum effieciency and minimizing decoherence. At Yale, 
Wallraff et al. coupled a charge qubit to a cavity made 
from a coplanar waveguide (CPW)[6j. In this experi- 
ment, a cavity with a high quality factor (Q) was used to 
enhance the dispersive phase shift of detected photons. 
However, using a high-Q cavity severely limits the read- 
out speed, naturally leading us to consider using a low-Q 
resonator instead. That said, low-Q reduces the disper- 
sive phase shift and implies the need for large detuning 
between the cavity and qubit, because qubit relaxation is 
strongly enhanced if the qubit frequency lies within the 
cavity resonance. Bertet et ai.jjj coupled a flux qubit to 
a low-frequency SQUID oscillator. Unfortunately, they 
found that the thermal photon noise of the oscillator, 



present even when not measuring, limited the qubit co- 
herence. 

We attempt to harmonize these apparently contradic- 
tory requirements. The standard methods of cavity QED 
[lj,[lj| are not appropriate for this analysis because they 
only treat the high-Q limit. Therefore, we take a new ap- 
proach, based on the quantum network theory introduced 
by Yurke and Denker |13| , which gives analytic results for 
arbitrary values of Q. Indeed, we find that the read-out 
is quantum limited (77 = 1) independent of Q. We go on 
to calculate the thermal dephasing time in the absence of 
measurement, t° . We find that low Q can compensate 
the effects of thermal photon noise on i° by suppress- 
ing low-frequency photon fluctuations. We also find that 
large detuning compensates for a the smaller phase shift 
by accommodating stronger driving. We go on to discuss 
optimal parameters for a realistic single-shot read-out. 

Our approach is general, and we use it to derive both 
the Hamiltonians for the measurement of the effective 
capacitance (inductance) of a charge (flux) qubit. How- 
ever, we focus the optimization on the Cooper-pair box 
operated as a charge qubit 01, capacitively coupled to 
a low-frequency, lumped-circuit LC-resonator. In this 
setup, the effective capacitance of the Cooper-pair box 
was recently measured by two different groups 
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FIG. 1: A CPB (flux) qubit capacitively (inductively) coupled 
to an LC-oscillator, which is attached to a transmission line. 

The Cooper-pair box qubit consists of a small super- 
conducting island coupled to a superconducting reservoir 
through a Josephson junction characterized by its capac- 
itance, Cj, and Josephson energy, Ej. The island is 
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coupled to a gate voltage, V g , through the gate capaci- 
tance, C g , and to the LC-oscillator via the coupling ca- 
pacitance, C m . The oscillator is in turn coupled, via C c , 
to a transmission line, here modelled as a semi-infinite 
line of LC-oscillators. The transmission line is charac- 
terized by its capacitance, C*r, and inductance, Lt, per 
unit length. (See Fig. [H) 

Since we have a Josephson junction in our circuit, we 
choose as our coordinates the phase differences across the 
circuit elements <\> a it) = j dt'V a (t'). In the transmission 
line, we number the phases across the capacitances start- 
ing from the oscillator, giving the classical Laerangianflij 

L = ^^-^c4>j-C g V g h + EjcoBcl>j + 
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where C q b = Cj+C g +C m , C osc = C+C m . Conjugate to 
each coordinate <p a , the momentum q a has the dimension 
of charge, and [<t>a., <lp] = ihS a f3- 

During read-out the oscillator and transmission line 
charges {qc and qi) will be driven harmonically at a 
frequency close to the bare oscillator frequency ojq = 
1 / ^jL(C osc + C c ), which we consider much lower than 
the qubit frequency w q b = Ej/h. The amplitude of the 
charge oscillations induced on the qubit island is deter- 
mined by the driving strength, which we characterize by 
a parameter (3 through 



(qc(t) + qi(t))r 



= /3 



Ej 



(2) 



where Ec — e 2 /2C q b is the charging energy of the qubit 
and k — C m /C osc <C 1. For f3 < 1, we can neglect 
Landau-Zener transitions and the qubit will follow the 
oscillator dynamics adiabatically. 

At the charge degeneracy point, we rotate the full 
Hamiltonian to the qubit eigenbasis (a x «-> <r z ) and ex- 
pand the eigenenergies to second order in the induced 
charge fluctuations 
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where gc = %n 2 E 2 1 /e 2 Ej indicates the qubit-oscillator 
coupling. From the qubit perspective, the slow trans- 
verse charge fluctuations induced by the oscillator re- 
sult in second order longitudinal fluctuations as described 
by the (qc + qi ) 2 cr z -term 0. The same term causes a 
state dependent shift in the oscillator's electrostatic en- 
ergy. This may be modelled as an effective oscillator ca- 
C osc ± gcC 2 sc = C osc ± C Q , where 
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the last term is the state dependent quantum capaci- 
tance Cq. This in turn gives different oscillator reso- 
nance frequencies for the qubit in the ground/excited 
state ujf/ c = lo (1 =p Cq/2(C osc + Cc)), which was re- 
cently experimentally measured in Refs. |lllll2j. 

The same type of analysis can be done for the flux 
qubit in Fig. ^16]. The oscillator part of this Hamilto- 
nian is 
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where now the qubit-oscillator coupling is <?l = 
2M 2 (4>j) 2 /AL 2 L 2 . Here (<j)j) is the average flux gen- 
erated by the circulating current in the Josephson loop, 
M is the mutual inductance between the oscillator and 
qubit, L q is the qubit inductance and A is the qubit 
level splitting. In this environment, the Josephson junc- 
tion acts as an effective inductance L 9 ic = L ± L 2 gL 
giving the state dependent resonance frequencies ui 9 J e — 
LO {lTLg L /2). 

The equations for the transmission line coordinates 
and momenta, in the limit Ax — > 0, correspond to 
the massless Klein- Gordon equation, which has traveling 
wave solutions, (f>(x,t) — 4>i n (t — x/v) + 4> ou t(t + x/v), 
with components moving towards and away from the 
oscillator with a velocity v = \/\JLtCt- The driving 
source determines the in- field, 4>i n {t — x/v), from which 
we may derive the detectable out-field, 4> ut(t + x/v), as 
well as the charge field q(t) = qc(t) + qi{t) and phase 
field 4>c(t), giving the backaction on the charge and flux 
qubit, respectively. In so doing, we solve the Heisenberg 
equations of motion corresponding to the Hamiltonian in 
Eq. J^J. Since they are linear, this is straightforward in 
the Fourier representation. The linearity is a consequence 
of the weak coupling to the qubit, which is the only non- 
linear circuit element. However, we are free to consider 
any strength of the oscillator-transmission line-coupling. 

First analyzing the charge qubit read-out we de- 
fine the state-dependent transconductance x s ^ e ( w ) = 
q 9 / e (lu) / 4>i n {ui) , and write the general solution 

o/ P , \ i2C c Cosc Lu 3 

x 9 ' M = 



1 - (C 9 t + C c )Llj 2 - iuC c Z (l - L^Cl'sl) 



where Zq = ^JLt/Ct is the characteristic impedance of 
the transmission line. Since there is no dissipation in the 
lumped circuit, we have = |<^ n (cj)|, and the 

reflected signal is characterized by a frequency depen- 
dent phase shift Lp 9 J e (uj). By introducing creation and 
annihilation operators for the in-field satisfying the com- 
mutation relations [a u , a^,] = 8{u> — u>') and [a u , a u >] = 0, 
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we can write 



<t> in {x,t) = N jH ' ^(ye-Mt-xh) + h .c), 
q 9/e (t) = N J~ '^L (x 9,e ^)a u e-^ + h.c) , (6) 



with normalization N — \JUZqJ4k. We model a sig- 
nal generator by setting the in-field to a coherent state 
with a narrow-band distribution around the drive fre- 
quency LUd and an intensity corresponding to r m pho- 
tons per second. Using the standard model of homo- 
dyne detection, the time t ms needed to distinguish the 
phase difference 6ip r = ipt(tVd) — ^(wd) in the out-field 
is t~g = 4r in sin 2 (Sip r /2) [15j. We can then compare 
this to the measurement-induced qubit dephasing. In 
the weak qubit-oscillator coupling regime, we straight- 
forwardly calculate the decay of the off-diagonal ele- 
ment of the reduced qubit density matrix |((7 + (i))| = 
\(a + (0))\e~ t / t «> and find the dephasing rate 
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The quantum efficiency r/ can now be evaluated for ar- 
bitrary circuit parameters. First we note that the drive 
strength Ti n cancels in the expression for r\. The LC- 
oscillator has a quality factor Q = (C osc + C c )/C^Z ujo, 
and close to resonance we can approximate Eq. © as 
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Driving the circuit at the average resonance frequency 
LUd = (wg + Wq)/2 the quantum efficiency is 
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where x = QCq / (C osc + C c ). Somewhat surprisingly, we 
find a quantum efficiency of unity independent of Q and 
the coupling to the qubit Cq /C OS c- For non-optimal drive 
frequencies the efficiency is below one. The inefficiency is 
related to the information stored in the state-dependent 
delay time of the measurement pulse. 

For flux qubit read-out Eq. J7J apply with the substi- 
tution gc — > <?l and using the relevant transfer function 
x(w) 9 / e = (f>£, e (uj)/(j) m (uj). In this case we also find full 
quantum efficiency independent of Q. 

We can now relate this result to the quantum effi- 
ciency of some other read-out methods. Detecting the 
phase of the signal transmitted through the CPW cav- 
ity in the Yale experiment has half the optimal efficiency 
r\ = 1/20|, where the inefficiency is related to the re- 
flected signal, which is not detected. In our system, the 



efficiency is optimal since the whole signal is reflected 
and detected. In the point contact detector, electrons 
in a range of energies, as specified by the driving volt- 
age, are used to probe the contact. Thus a particular 
energy-dependence of the point contact transmission is 
needed for the efficiency to approach unity I 11 our 
system, the signal source is essentially monochromatic, 
probing the system at a single energy. Thus, there are 
no constraints on the energy dependence of the scatter- 
ing, allowing, e.g. for a \ow-Q oscillator. 

Finding full quantum efficiency for a wide range of pa- 
rameters leads us to consider optimizing the circuit pa- 
rameters for other criteria. One important figure of merit 
is the rate of qubit dephasing induced by the measure- 
ment device when not measuring. In this case, we assume 
a thermal distribution with temperature T of the in-field 
and find the dephasing time 
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where n(wo) = \/{e hw °/ kBT — 1) is the thermal photon 
occupation number at the oscillator frequency. This ex- 
pression is similar to what is found for a superconducting 
flux qubit coupled to a DC SQUID 9]. To minimize off- 
state dephasing one should cool the oscillator, but it is 
also advantageous to reduce the oscillator-qubit coupling 
and to use a low-<5 oscillator. Basically, the integral of 
the spectral density of the fluctuations, the variance, is 
fixed by thermodynamics, independent of Q. Therefore, 
lowering Q stretches the same power over a wider band- 
width, reducing the low-frequency power. 

Another figure of merit is the relaxation time of the 
qubit, fx, induced by the measurement device. We evalu- 
ate t\ using standard weak coupling expressions 0, giving 
that it is proportional to the real part of the impedance 
seen from the qubit, at the qubit frequency uj q b- For 
the simple oscillator circuit shown in Fig. ^ with a low 
Q, t\ could be as short as a few hundred nanoseconds. 
The solution to this apparent problem is to add a non- 
dissipative low-pass filter between the transmission line 
and the oscillator. Due to the large frequency difference 
between the oscillator and qubit, this is straightforward. 
In principle, commercially available pi-filter can improve 
t\ by a factor of 1000 while not affecting the low fre- 
quency measurement properties. This is well beyond the 
point where other sources of relaxation will dominate. 

Along with these intrinsic figures of merit, we should 
also consider the performance of the read-out constrained 
by existing technology. In particular, we want to evalu- 
ate the measurement time given that we first amplify 
the reflected signal with an amplifier characterized by an 
impedance, Zq, and a noise temperature, TV, which ex- 
ceeds the quantum limit. We also assume the input of 
the amplifier is coupled to a matched load, e.g. using 
an isolator. Now, the maximum number of photons in 
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FIG. 2: The number of oscillator photons (n max ), off-state 
thermal dephasing time (£^), oscillator quality factor (Q), 
and phase shift (ip r ) as a function of read-out time. 

the oscillator, n max = C osc (ef3Ej /2E c C m ) 2 /Tiuj d , is de- 
termined by how hard we can drive the qubit, which is 
limited by the width of the quantum capacitance peak in 
gate charge. Roughly, we can use the half width at half 
maximum (HWHM) of the peak, giving f3 = \/2 2 / 3 — 1/2 
in Eq. J5J. Taking the phase shift between the ground 
and excited state as (p r — arctan(4<5CQ/C osc ), we can 
estimate the measurement time as 

_„ n amp Q Z kBT]yCc 

where n amp = ksT^ /Tiujd is the number of noise photons 
in the amplifier, A = C m /C q b, and we have assumed a 
signal-to-noise ratio in phase of 2, which corresponds to 
the same definition used for the quantum limit. In the 
last approximation, we assume ip r ~ ^QCq/C osc . 

We then want to find a circuit design that gives the 
fastest measurement time while still protecting the qubit 
from thermal dephasing. In addition, we require that 
the measurement time be longer than the response time 
of the resonant circuit, tq = Q/luq. This ensures that the 
reflected measurement pulse is not distorted and that its 
delay is not state dependent. In Fig. 2, we plot opti- 
mized values for n max , <p r , Q, and as a function of 
t ms - We assume t ms /rg = 4, T = 20 mK, A = 0.5, 
w /2tt = 650 MHz, u) qb /2w = 5 GHz, Z Q = 50O and 
Ttv = 0.8 K which implies n amp w 25. We see that 
measurement times of order 50 ns are readily achievable 
while still maintaining off-state dephasing times greater 
than 10 fis. 

We can compare this result with the estimates of ^ij , 
trying to understand how the estimates here are 2-3 or- 
ders of magnitude faster. Fundamentally, the speed-up 
comes from having a much lower Q, which determines 
the response time of the resonator. However, lowering 
Q alone also decreases ip r . In fact, if nothing is done 
to compensate this, the measurement time actually in- 
creases as 1/Q. The way to compensate is to increase 



nmax- This is where the large detuning of the our pro- 
posal is key. In the near-resonant design of |l4| . n max is 
limited by the fact that the shift of the qubit frequency 
due to photon occupation must be less than the detuning. 
If this condition is violated, qubit relaxation is strongly 
enhanced by the cavity. This implies n max » 100 and 
n a mp/n max ^>r ~ 1- This is not a limitation for a low- 
frequency oscillator, however, allowing for much harder 
driving. For the the case of t ms — 50 ns above, we have 
in fact n osc s» 1200, maintaining n amp /n max ip 2 w 1/2 de- 
spite the small phase shift. In fact, this ratio is constant 
for the whole range of plotted parameters. 

In summary, we have discussed how to optimize a gen- 
eral dispersive qubit read-out for speed, while maintain- 
ing maximum quantum efficiency and long dephasing 
times. We have found that reading out a charge qubit 
by measuring its quantum capacitance, through the fre- 
quency shift of a lumped circuit LC-oscillator, is quantum 
efficient (rj = 1) independent of the oscillator quality fac- 
tor. We also find that a low Q protects the qubit from 
thermal dephasing in the off-state, while increasing the 
speed of the read-out. Furthermore, we have discussed 
how to optimize the circuit parameters to achieve single- 
shot read-out using a commercial amplifier. 
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